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ABSTRACT 


Global temperatures are obtained using a linear least 
Squares regression method with satellite radiation measure- 
ments, particularly those obtained from the SIRS-B aboard 
the NIMBUS IV satellite. Regression equations relating 
temperature to spectral radiance observations are employed. 
The regression equations were then applied to independent 
observations to determine the feasibility of temperature 


determinations over Sparse data regions. 
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L. INTRODUCTION 


The launch of the NIMBUS IV on 8 April 1970 was the 
second step in a process of developing a means of determin- 
ing the three-dimensional state of the atmosphere by indir- 
ect means. The orbit of the NIMBUS IV was a near-circular, 
polar-orbit, with a height of approximately 600 nm. The 
new instrument, Known as SIRS-B (which is similar to the 
instrument on NIMBUS III, known as SIRS-A)) has the added 
capabilities of "side looking" scan for greater area cover- 
age, and the determination of the tropospheric water vapor 
Eentent. 

The calibration procedure used to transform the radia- 
tion sensed in the eight channels of the SIRS-B to spectral 
Padiances is discussed in the NIMBUS IV USERS GUIDE (1970). 
The temperature profile of the atmosphere within the field 
of view, comprising an area of ianye at the subsatellite 
Petme, is deduced using the channel spectral-radiance en- 
semble. The particular wavenumbers sensed are small inter- 
vals in the atmospheric window and the 15-micron band of 


co,. The channel numbers and wavenumbers are listed: 


CH. NO. 1 2 3 4 3: 6 7 8 


WAVE. NO. 899.0 750.0 734.0 709.0 701.0 692.0 679.8 666.7 
v, Cem™?) 








Through the equation of radiative transfer, the meas- 
ured radiances are related to the temperature-structure 


T(p) of the atmosphere, the temperature at the earth's sur- 


face and that at the cloud tops which may be within the 
field-of-view. Assuming that a single cloud layer exists, 
the following set of integral equations can be constructed 


for the measured radiance Nv.) inmeehamnel i: 


Cc 


N(v,)=A{BIv,, TOP] TO oP.) - Sf “Bly, .T(p) ]dt (yy sp) | 


. (1) 
S 


+ (1-a) { BIv,T(p_)] t(v,,p.) -,f Blv,,T(p)l]dt(v,,p) | 


iP so. wee so Loethe channel number. 


N(v5) spectral radiance at wavenumber V5: 


Planck radiance at wavenumber v, and tempera- 


BLv.,T(p) ] i 


ture, (Eq. (3) below). 


T(V,5P)* fractional transmittance of the atmosphere 
in the spectral interval centered at wave- 
number Vs and assumed known as a function of 
the pressure p at the radiating level. 

Me mene product of sthe fraction of cloud cover 


within the field-of-view and the cloud emis- 
sivity (which is usually assumed to be unity). 
Subscripts c and s refer to cloud-top and earth surface. 
Meesemopecial case of equation (1)) where A=0 gives the 
clear-column or corrected radiance NOv3)- This term is 
contained within the second brace of equation (1). It is 


this quantity which must be determined in all channels in 
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Demeerbowconvert the corrected radiance set into a vertical 
temperature profile T(p). The clear-column spectral radi- 


ance at the top of the atmosphere at wavenumber v, is there- 


fore: 
Nv) = BIv,,T@p,) ITO, >P.) 
*s (2) 
= id BIv,,T(p) ] dt(v.,P.) 
where 
3 , 74 
Ei ee eileen, exp (cv, /T(p)} = 1} (3) 


is the Planck function at wavenumber Va 
Here T=T(p) is the temperature (assumed to be in thermo- 


dynamic equilibrium at level p), and c, and c, are known 


1 Z 
constants 


1.GOGl oe erg ane gee" Coa 


= 


1.43868 cm °K 


2 
TCV. »P) is the transmittance at wavenumber A of 
the atmosphere above level p. 
Smith, Woolf and Jacobs (1970) have used equations (2) 
and (3) to define a weighted-mean black-body temperature 
te a in terms of the corrected radiance Nv ,) which would 


be observed from an atmosphere of structure T=T(p) 


N(v,) 3 Bee Da) (4) 
Here Tae is independent of pressure. Solving the relation- 
sauaps (3) and (4) for Tai gives 

” 3 
mG) = € v,/in{ fe,v, /N.(v,)1 +1} (5) 


Smith et al (1970) have provided multiple regression equa- 


tions which relate NMC effective grid point temperature 
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wai at SIRS-A scan spots to the sixteen predictors T, OY) 
and on). They have done this by relating a two-week 
NMC data file of temperatures at standard pressure levels 
P, to the corrected radiances, and arrived at an equation 


of the form 


8 
Tp s) = Tep5) +E atvyepy) ITO) He 
: = 2 
fE ahve s) [TgQyg~ TRO] (6) 


A similar study will be made for T(500) sea T(300) using 
a data-base of corrected radiances from SIRS-B and NESC de- 
rived temperature profiles T(P 5) together with simultaneous 
FNWC effective grid-point soundings (available from the 
Fleet Numerical Weather Central ADP tapes) at locations of 
the scan-spots in both space and time. For the SIRS-B data 
file, a revised NESC formula similar to equation (6) was 
available for the conversion of corrected radiance-sets 


into a statistically consistent profile T=T(p). 


a7 





Il. DATA PROCESSING 


A. THE ORIGINAL DATA 

The original data utilized in this paper were provided 
by the National Environmental Satellite Center (NESC), and 
consisted of computer listings of a selected set of SIRS-B 
weeecopoets recorded during December 1970. in addition to 
the geographical coordinates and time of each observation, 
there were also listed the raw (uncorrected) radiances in 
each channel. The raw radiances are assumed to be derived 
from equation (1) for the one cloud-layer model and the 
corrected radiances resulted from the systematic elimina- 
jmion of cloud epee na naked effects from each channel [fol- 
Heowineg a technique of Smith et al (1970)]. 

The atmospheric window channel (channel 1) was inopera- 
tive during this period; therefore, the values of the cor- 
rected window channel were chosen from locations where the 
Surface temperature, based on the NMC data file, was known 
to be at a pressure-level close to p = 1000 mb, where T is 
known. Furthermore, since the radiance-correction computa- 
tion involves iterative use of equations (2) and (6) (using 
an equivalent "clear-column" temperature profile, T=T(p), 
at each iteration), all original radiance-sets chosen could 
thus be "tied down" by the selection of scan-spots having 


surface pressure values Pa LOO 0 mb with tT. known. 


= 


13 





It should be noted that the data samples provided were 
enhanced by the side-viewing capability of the SIRS-B. The 
side-viewing capability enabled the scanning of more atmos- 
phere-earth columns both on and off the subtrack for which 
the approximation ea 1000 mb was valid. Slant-columns 
thus viewed are useful for reduction to clear-column tem- 
perature profiles T = T(p) after the slant path is normalized 
to a single atmospheric depth over the spot viewed. 

In general, virtually all scan spots have some effective 
cloud cover (A > 0) in view of the large area (225 km) be- 
amg viewed. Of the seven CO, channels being used, the two 
channels subject to the maximum attenuation by the inter- 
mosttron of a “cool?” cloud top within the desired T = T(p) 
profile are channels 2 and 3. 

Smith, Woolf, and Jacobs (1970) have indicated the pro- 
cedure involved for correcting for cloud-contamination of 
the eight channels of radiances. Based upon the actual NMC 
climatological update of each of the standard pressure- 
Mergeno (10, 30, 50, 100, 150, 200, 250, 300, 400, 500, 

700, 850, and 1000 mb) they have derived multiple regression 
equations relating ee to the clear-column-radiance 
N(v,), or to the blackbody temperature-equivalent of the 
clear-column radiance. For example, Smith et al (1970) 
obtained regression equations essentially of the form of 
equation (6) and have listed the least-squares best-fit 
coefficients obtained for a period in September 1969 appli- 


cable in a zonal band 35-55°N. As a second approximation, 
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Smith et al have also generated a six-predictor analog of 
equation (6) with channels 2 and 3 deleted for the same 
time frame as utilized for the eight-predictor equation. 
In both cases, they have deduced equation (6) and its six- 
predictor analog by employing a large sample-size of simul- 
taneous values of corrected radiance-sets and of values of 
ar) (at standard pressure levels), interpolated from the 
WHC £Lile to the location of the SIRS scan spot. The six- 
predictor and eight-predictor regression equations are con- 
tinually updated to a revised two-week data base after 
three additional data-days are accumulated. Such equations 
were also available during the SIRS-B period of December 
70. 

mite precedure £or correction of the ith channel radi- 
ance makes iterative use of the updated SIRS-B operational 


six-predictor regression equation in deducing a first guess 


of T=T(p). At the first guess, the uppermost five radiances 
(Nis eeeoNg) which are least affected by cloud contamination 
are used as the first guess correct-radiance set, together 


with the prescribed N, from the window-channel blackbody 


i 


' at the scan spot. With these six-predictors known, a 


Pfs,’ 
first guess of the temperature profile, ae at all 
standard levels is computed. After the first-guess of the 
temperature profile is made, the profile is then applied to 
channels 2 and 3 to determine the "best fit" for cloud 


height and effective cloud cover A using equation (1), re- 


written for simplicity as 


Mes 





MT) PAY N04 T@) pl) + AY [v,,1(p),P) (7) 
Here Y represents the contents of the first brace of equa- 
tiom (1). Also c(v,) is defined as 

c(v,) = Ne (v,) = N(v,) (8) 
for the correction to be applied to the radiances for chan- 
nels 4 through 8. Writing equation (7) for channels 2 and 
S7tnterms of the corrected vahues of No based on T(p) (the 
clear-column temperature) and the corresponding computation 
of Y from the top of the clouds gives two simultaneous 
equations. From these two simultaneous equations, one 
solves for the two unknowns, A and Po» USiNe Bindi. eet 
methods for solving for P. (i.e. making separate tests for 
a ZOO 2 SO. 5007) 850 mb) and, in turn, determining the 
value of A which "best-fits" equation (7) expressed in 
terms of N(v,) and N(v,). These values of A and Pp, are 
Substituted into equation (7) for a eacond corrected esti- 


mate of N ao: eee 


4? and a new temperature profile is calcu- 


8 
lated. The iterative procedure is continued in this manner 
and usually converges within four iterations, i.e., the 
correction c(v;) of (8) usually converges 

ees CCU Li mat ) 

i at 

Honma SuUlLtable tolerance Within four iterations. At Ehis 
point, the finalized temperature profile T(p) may be used 
to compute the clear-column radiances for channels 2 and 3 


as well as the five uppermost CO, channels using Eq. (1). 


Z 
Cloud influence is thus compensated when the final choice 


of temperature profile has been made. 
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In a sizeable fraction of the cases selected, the con- 
vergence c(V,) required more than four iterations which 
means that a relaxed requirement for convergence hel to be 
employed. This generally occurred with cloud tops at or 
above 500 mb. Such cases of below-normal attenuation re- 
covery were indicated on the NESC data-listings as decreased 
reliability. Nevertheless, no data cases were discarded 


for such reasons. 


B. ADJOINING OF FNWC DATA 

Analyzed temperature-field data for 500 and 300 mb cover- 
ing the period from 00 GMT 24 December 1970 through OO GMT 
29 December 1970 was provided by Fleet Numerical Weather 
Central (FNWC). The FNWC temperature data had not been 
correlated with radiances as was the NESC data. The FNWC 
temperature fields at 500 and 300 mb (T.(5) and T,(3)] were 
spatially and temporally Lae bILARed Be obtain consistency 
Mite the SIRS soundings. Spatial interpolation was accon- 
plished using Bessel's central-difference interpolation 
[Maltiner, 1971] to the latitude and longitude of the SIRS-B 
soundings. Temporal interpolation to the time of the scan- 
Smet Was accomplished utilizing linear interpolation between 
temperature fields preceding and succeeding the time of the 
SIRS sounding. These methods were employed to obtain 
statistical dependence of T, 05) SLi T, G3) upon the corrected 


radiance set N are 


yore This necessitated the printout of 


8° 


ye 





FNWC interpolated temperatures in a format consistent with 
the requirements of the program BIMED O2R, to be used in 


CTHer ) . 


C. DATA DESIGNATION 

The original NESC data comprised some 207 to 329 scan- 
Spots per day and were identified by the specific time of 
the SIRS soundings and their latitude and longitude. The 
Original FNWC data tapes consisted of temperatures at the 
lewels 300 and 500 mb, at each grid point of the 63 x 63 
Northern Hemisphere FNWC grid. The FNWC data at scan spots 
were obtained by interpolation as previously described. 

A three-day data base of both NESC and FNWC data was 
thencompiled using data from 0600 GMT to 1800 GMT for 24, 
25, and 26 December 1970. <A four-day data base was com- 
piled by the addition of data from 0600 to 1800 GMT 27 
December 1970 to the three-day data base. Data sets for 
independent regression tests were Sonoiled using data cen- 
tered on 00 and 12 GMT 27 December and OO and 12 GMT, 28 
December (i.e., data-base times plus 12 and 24 hours). 

All data sets were divided into three bands. These 
bands were bounded by latitudes: 20-40°N, band 1; 35-55°N, 
bamd 2: and 50-70°N, band 3. The bands were overlapped to 
reduce the possibilities of discontinuities across the 


boundaries when regression calculations were made. 
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The designation of dependent and independent variables 
was made and the dependent variable was designated T(i,j) 


and the independent variables were designated as shown: 


on. HO. 1 Z 3 4 5) 6 7 8 
wee ee onmeeo 7. 0750.0 734.0 70920 701.0 692.0 679.8 668.7 
VAR. NO. Ny N, Ny Ny N. Ne No Ne 


For each scan spot selected in the test period, the 
time and location were encoded on IBM cards, as were the 
temperatures T05) 5 T(3); TL 5); T, (3). Here, the sub- 
script N indicates an NESC-generated temperature, using the 
convergent procedure described previously, and the subscript 
F denotes interpolated FNWC temperatures. Also encoded 
onto IBM cards Pee a eh @lghe channels of corrected radi-— 
ances(provided by NESC using the method described in II.A.) 
for each of the scan spots. All data were encoded in a 
Manner consistent with the requirements of the stepwise 
regression program, BIMED O02R (Dixon 1966) which will be 


G@usccussed later. 
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ite tie REGREoos LON Merl HOD 


me We SLEPWISE REGRESSION PROCEDURE 

Many problems in research require the extensive analysis 
of large amounts of data. The data handling process should 
be" mate “as "autonatic~and~rapidwasrwpossible. »«Bhe «appropriate 
tools for analyzing large data samples with numerous inde- 
pendent variables often require multivariate statistical re- 
gression techniques to attack the problem. The Biomedical 
Computer Programs (Dixon 1966) were developed to provide 
the techniques necessary for the statistical and mathemati- 
cal analyses required. One of these programs is the Step- 
wise Regression Analysis Program BIMED O2R. 

BIMED O2R computes, in a stepwise manner, a sequence of 
multiple linear regression equations. One variable is added 
to the regression equation at each step. The variable added 
is that which accomplishes the greatest reduction in the 
previously unexplained sum of squares. This variable is 
also the one which has the highest partial correlation with 
the dependent variable at the particular step in the analy- 
sis of the variance. Equivalently it is the variable which 
would have the highest F-value if it were added. The F- 
value, or F-statistic upon entry, Fy is expressed at step k 


as (Dixon 1966) 


TCsE Velie LOGE ak 1) 
TCE) 


FY (1,n-k-1) 
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where 
&(C.E.V.k) is the percent cumulative explained variance, 
step k. 
&4(C.E.V.k-1) is the percent cumulative explained 
variance, step k-l. 
&(U.E.V.k) is the percent unexplained variance remain- 
ire at step k. 
This study utilized a statistical model expressible in 
the form 


T = © TON TOONS TC QNA TC LN, TCLN TC N FC IN FCN, (9) 


where oR through C, are to be determined by the stepwise 


8 


least squares technique. Equation (9) may be expressed in 


tire more compact form utilizing the matrix-product equation 


Na 
ee eae) (10) 


In (10), the notation T. (355) will be used when the NESC 
data at 300 or 500 mb (respectively) are correlated against 
the predictors (N,> rae »NQ)3 whereas T, (3,5) is the nota- 
tion for the corresponding predictands with FNWC tempera- 
tures. Equation (10) may also be used in application with 
independent data occurring up to 12 or 24 hours later than 
the dependent sample which yielded the particular coeffi- 


cient matrix G C reer 
( ond 1 > 


wesc bE STATISTICAL PARAMETERS 
Available as outputs from the BIMED O2R program are 
“several related statistical parameters. Among these sta- 


tistical parameters are: 


oy 





a. multiple R 

Eeeestandard Crror of estimate, S./E: 
c¢. mean value, ig 

d. standard deviation o 

e. F-value 

ts R? 


Mme significance Or R? in@rékrarrvonsiiap™to-S ..E. is as 


follows (Crow et al 1955): 


eate@ie 10; (eke) Gene (11) 


(suzy)? 
where 


as =e? 
o* =f Co 1) e/a 


is the variance of the temperature sample, and 


n = sample size 
i = sample-element index 
k = number of predictors 


Purethermore, R is the multiple correlation coefficient after 

k predictors are added, with emphasis on k=8, in order to 

derive maximum information from the SIRS-sounding. It fol- 

lows that the percentage unexplained variance ers) f£ Ons 

sample sizes in the range 171 to 228 as contained in the 

dependent data sample are very closely approximated by 
(omeo)- = ee 

The percentage explained variance is then given by ne that 

is 

ee (se Bo)" (12) 


since (n-k-1)/(n-1) is close to unity for the dependent 


oe 





sample sizes considered here. The factor (n-k-1)/(n-1) 
becomes (n-2)/(n-1) = 1 for the independent test cases 
considered since in this case, there is only one way of 
forming the predictor to be listed. 
ett eekeGRESSLON PROCEDURE APPLIEED TO THE DEPENDENT 
SAMPLES 
The test on the dependent data was made using the least- 
Squares, stepwise regression technique of BMD O2R to develop 


the coefficients Coe + Ce. The coefficients for the 


rr? 
three-day and four-day dependent samples for the T(5) and 
T(3) regression equations were computed by use of this 
technique. This test was performed on both the NESC-regres-— 
sion generated data and the interpolated FNWC data. The 
coefficients for the NESC three-day data sample are listed 
femable 2, part (a) for 500 mb; and part (b) for 300 mb. 
Table 2 contains similar results for the four-day sample 
base of NESC data. The coefficients for the FNWC data are 
listed in Table 3 (three-day case) and Table 4(four-day). 
The relevant statistics for the tests on the several 
dependent data bases are listed in Table 5 and 6. The 
statistics included in these tables are the mean, standard 


; Z ; 
deviation, R ', and standard error. A discussion of the 


Sevtistics of prime interest is provided in IV. 
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TABLE 1. Regression coefficients for NESC dependent data. 


(a) 500 mb, three-day sample 

CONST 20-40N 35-55N 50-70N 

Co -43.1058 ee -90.9234 
Ci -0.0327 O44 37 C215 08 
C. =-0.0932 -0.4795 -0.2562 
C. G2 s922 -0.2670 0.0 

C), 0.6949 S002 52 259128 
C. -1.8871 -2.6461 — ees 
Ce -0.9188 =)07 26 -0.1134 
Cl 0.4129 OO. 1608 C2235 0 
Cy 0.6664 0.0764 -0.1635 

(b) 300 mb, three-day sample 

Cy = Or. 3 -113.6468 -108.4219 
Ci 0.4109 0. 5032 0.3409 
Co eO.29 55 -0.3192 -0.2337 
C. -1.4766 -1.4030 -0.6563 
Cy 3.0236 5211s WE Seats: 
C. -0.7404 OR OG57 nO 933 
i. -0.5295 0.07 0.1261 
Co =O 5055 SJ. 5 Oe 2S) -1,.1938 
C. 0.6865 0.4809 Os 922 


indicates that this variable (*) was 


not accepted 


because it failed to explain appreciably 
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TREBLE 2. 


Regression coetiicients for NESC dépendent data. 


(a) 500 mb, 


four-day sample 


CONST 20-40N 35-55N 50-70N 
Co -48.6906 = Ole oZ oe -92.4865 
Cy 0.0 0.4842 0.1143 
C., -0.0907 -0.6454 -0.1644 
C. 0.8032 -0.1472 -0.0980 
Cy, 0.6590 3-200 5, 0117 
C. -1.5652 =—2.6431 -1.7844 
o -0.9017 0.0" -0.1327 
CL 0.2990 Om 60S 0.1894 
C. 0.6797 0.0354 =O oa 
(b) 300 mb, four-day sample 
Co -101.7840 -108.5859 -108.1667 
Cc, Oneeze Ss U.5599 ONS 57 
C. -0.0794 -0.5138 =-0.3326 
C. -0.2788 -1.0924 -0.5185 
Cy, Ons 2ro ID Jf), 2) 2 Te AS 107) 
C. £22700 0.1432 1646 
C. -1.0121 O08 O1026 
Co -0.5885 -1.0406 =I 2326 
C. e295 Oy LS 0.4386 


* 
indicates that this variable 


(*) was not accepted 


because it failed to explain appreciably 
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TABLE 3. Regression coefficients for FNWC dependent data. 


(a) 500 mb, three-day sample 


CONST 20-40N 35-55N 50-70N 
7 -15.6121 -55.3823 -~85.1272 
S 0.3605 0.4119 9.0638 
., Sones. -0.8714 Oe 72S 
C. 0.7760 0.9170 Torey 3 
C, 1.2602 1.6287 0.6259 
G. -1.8509 -~1.9698 -0.6344 
CP =O ce os 0.1453 -~0.1520 
Ca 0.2310 -0.1858 0.0% 

e 0.4440 0.1834 -0.0545 


(b) 300 mb, three-day sample 


Co -56.7715 -107.1405 -99.6836 
ea 0.4459 0.2963 0.3651 
Sa -0.6807 0.0459 -0.6303 
ce -0.3890 -1.1144 0.3020 
c, 1.9185 2.1166 0.5930 
ce -0.5798 0.3357 0.9944 
Ce -0.9211 -0.1590 0.0444 
C. -0.2675 -0.5875 -0.6276 
Ce 0.5916 0.2311 0.0392 


* e 
indicates that this variable (*) was not accepted 


because it failed to explain appreciably 
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TABLE 4. Regression coefficients for FNWC dependent data. 


(a) 500 mb, four-day sample 


CONST 20-40N 35-55N 50-70N 
z -15.3487 -~51.3814 Scemeos 
om 0.4477 Oaso 72 0.0260 
Gis -1.0928 See =i) Se 
om 1.0306 1.3415 1.4995 
on plage 1.3383 0.8176 
ue RYT -1.9340 =) 5 ONE 
om -0.8981 0.2096 =) 1 Ay, 
e 0.2645 =Onsa6 -0.0979 
Gis 0.3767 0.1141 EO? 


(b) 300 mb, four-day sample 


Co -53.5409 -103.0841 =) Soe a 
Cy 0.4850 0.3744 ONG sie 7 
Cc, -0.7624 -0.1827 -0.6721 
C. (0 5 fe EA NO, =O 6659 O37 22 
Cy DO 2 Za LoS 0.5005 
Ce -0.6480 Omi 0b 7 7 
Ce SOR SO eots 0.01/71 
Cl -0.1607 -0.5045 (0) 7) 
Ce 0.4292 0.2066 O07 57 


+ 
indicates that this variable (*) was not accepted 


because it failed to explain appreciably 
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The detailed results of the tests on dependent data are 
available in Tables 5 and 6. However, the essential por- 
tions needed for comparisons by the different tests using 
dependent data have been listed in Table 7 where the values 
of fractional explained variance as originally listed in 
Tables 5 and 6 are summarized. Also listed in Table 7 are 
the values of R? averaged over the data-base specification 
period (or mode) for each band as well as the values of Re 
averaged over the bands for a given data-base mode. Com- 
parisons between regressions were based on the fractional 
explained variance a rather than on the values of standard 
errors, as ee has been normalized with respect to the sam- 
ple variances, which vary between bands. 

The third value listed in the band columns of both parts 
(a) and (b) of Table 7 is the average over the data-base 
modes in each band. The values listed in the mode-average 
row and band-average column are the overall averages by the 
statistical regression. Comparisons of R? are made of the 
overall averages and of band averages, as well as data-base 
mode averages at both the 500-and 300-mb levels. Differ- 
ences in values of certain compared R? are termed "% shrink- 

" 


age in the explained variance arising from the specifica- 


tion tests under comparison. 


=> 
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ma) Pew ARTSON OF NESC SPECIFICATION STATISTICS 

Examination was first made of the left hand side of 
Table / for comparison of the NESC regression Pees at 
300 mb relative to those at 500 mb. The outstanding 
Fesult of this comparison reveals R (3) to be smaller than 
Ri (5) by sO. 1660, or 16.6% when averaged over bands 1, 2, 
anda 3. Upon comparison of individual band results (averaged 
over the three and four day specification modes) one finds 
the same order of shrinkages distributed over the bands, 
Namely, an average of 16.6%, but with a larger than normal 
shrinkage in band 2 where the value of RY relleori Vio 
O78 57. 

Upon use of the four-day estimation equation C0) iia 
the respective cases (left side of Table 7) there was a 
small shrinkage in explained variance (relative to the 
three-day results) at both the 300- and 500-mb levels. 
This result was not quite anticipated since Smith et al] 
found it advisable to use data-banks of two weeks for 
specification. Nevertheless, the small differences between 
the RO for the two prediction-base methods, whether at 300 
or 500 mb, were not large enough to be considered signifi- 


cant in the tests conducted here. 


B. COMPARISON OF FNWC SPECIFICATION STATISTICS 
Here our comparisons will mainly concern the right hand 
Side of Table 7. Again, predictor-mode differences seem to 


be marginal so that for the present they will be disregarded. 
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Comparison of Re (3) with Re (5) shows that band average 
shrinkage in the amount Ra (5) ~ Re (3) = 0.1480 occurs. A 
breakdown of this statistic by bands shows pees in the 
ames 3 On, 22.394, and 12,874 for bands 1.) 25 and) 3. 
respectively. 

Comparisons made between NESC and FNWC regression re- 
Sults (left side and right side of Table 7) show approxi- 
Mately 2% difference between Ro (5) = 10.67 37o-and R (5) 
0.8957 and less than 0.5% difference at the 300-mb level. 
Thus, the FNWC specifications at both levels explained very 
nearly as much variance as did the NESC results. There 
were some variations between bands, but generally, the 
same geographic pattern persisted throughout: namely, 
moderate to high oa in bands 1 and 3, with the lowest value 
in ne occurring in band 2. This was true of all FNWC speci- 
fications (300 and 500 mb) as well as of the NESC 300 mb 
Speerrication. Ihe NESC 500 mb specification gave rise to 
an increase in R? with an increase in latitude. 

In summary, both NESC and FNWC regression equations 
tested well in terms of explained variance at 500 mb. Both 
showed diminished specification at 300 mb to values of 
ee 0.73. Finally, in no case was there a significant 
change when testing the four-day data base results against 


those for the three-day data base. 
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Meee 5. 


Dependent-data temperature specification 


Statistics resulting from pooled samples of NESC data. 


(a) 
Sample 
Band Size 
20-40N 7h 
35-55N 191 
50-70N deep] 
(bd) 
20-40N 206 
35-55N 228 
50-70N 216 
(c) 
20-40N 7 i 
35-55N iUieal 
90-70N 7 
(d) 
20-40N 206 
35-55N 228 
50-70N ZALO 


500 mb, 
Sed, 
Mean Dew. 
=o U2 9 5 Oro 730 
=—2 O52) EO 47 OS 
Seite O69 67.0397 
500 mb, 
=15.4369 “6.7276 
Son LO23 "0.47376 
—34.21760) 7.0932 
300 mb, three-day 
~42.7895 4.6668 
50.4869 4.7427 
-54.6667 3.9610 
3007 mp; 
=42,/230 44.0392 
-50.3509 4.6875 
-54.7778 3.9610 
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three-day sample base 


R2 


ONE (ers) S83 


OR ales: 


CVs Sian 


four-day sample base 


Ore Oia 79 


0.9040 


02 335°8 


sample base 


Oe? P13 


0.6842 


OR 655 


four-day sample base 


0.6548 


Ores 71 


037970 


Sitar 


2h or 


Deeley 3) 1118, 
22 


1.8724 


Dwell SD 
Zr 02 


0 S10 


Po a4 
Di 2 


1.8874 


2.7609 
220057 


es te do toils 








TABLE 6. 


Band 


20-40N 


bol ON 


07 0n 


20-40N 


35-55N 


90-70N 


20-40N 


B55 N 


pO-70N 


20-40N 


35-3 om 


50-70N 


Dependent-data temperature specification 
Statistics resulting from pooled samples of FNWC data. 


(a) 


Sample 
size 


17k 
191 
177 
(b) 
206 
228 


216 


728 


216 


900 mb, 
SIL 
Mean Dien a. 
“15.9520 6.6156 
a oa ey A a ee 2 
so daee soo. 7 3364 


500 mb, four-day sample base 


ao 7 99 Go 
-28.0130 7.1574 
-34.2405 7.3470 


300 mb, three-day 


—42.2974 5.2983 
-50.8692 4.6688 
-54.9912 4.3548 


300 mb, four-day sample base 


-42.3044 5.1558 
=50', 7638 4.6365 
=o. 299° 24.342 


2 


three-day sample base 


2 


R 


0267 60 


0.8704 


0.8816 


OR o6 a 


0.8598 


0.8894 


Sample base 


0.7833 
0.6413 


0.7494 


0.7747 


0.6411 


0.7641 


Std. 
Error 


2239861 


2.7045 


Hd 8 Ter] 


2eG232 
22 OS 


2.4842 
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VY. ° INDEPENDENT DATA TEST RESULTS 


Verification of the dependent regression equations was 
accomplished using the regression coefficients (described 
in III) and the 12-and 24-hour independent data samples. 
The equations of form (10) developed from both the three-day 
and four-day samples were tested. The independent tests 
generate equations of the general form 


DC S53) “= By a3 B,1065,3) C13) 


~~ 


wnere By and By are regression constants and T represents 
the temperature-estimators from the dependent regression 
equation (10) with .coefficients appropriate to the level, 
the band and the mode-stratifications (Tables 1,2,3,4). 

The regression analysis for these tests is similar to the 
techniques applied in IV in that T(5,3) is compared by re- 
gression methods to the observed temperature at each scan- 
spot. The results of the independent tests are given in 
Tapes 8, 9, 10, and 1l. Table 12 is constructed in a man-— 
ner similar to that used to develop Table 7. The major 
difference between Tables 12 and 7 is that the 12- and 24- 
hour independent samples are statistically pooled in appli- 
cation of the independent regression test in order to obtain 
a larger sample size for each band, and also for ease of 


comparisons. The results of the pooled independent data 


verifications are presented in Table 12. 


= 
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Internal comparisons of the values of fractional ex- 
plained variance for the independent data yielded results 
quite similar to those obtained from the dependent data. 
The comparisons are stratified in two manners; by bands and 


by modes, as well as by levels. 


A. COMPARISONS OF NESC INDEPENDENT VERIFICATIONS 

From the left side of Table 12, it is easily observed 
that the largest overall average Ree 0705977 occurred in 
connection with the 500-mb verification of Ty 05). The cor- 
responding mean NESC result at 300 mb verified with a per- 
centage explained variance R (3) =" 04.7 933 ieee oan o Wietene Fa 
shrinkage of 8.59% relative to Ri (5). 

As noted in IV(A and B), both NESC dependent specifica- 
tions (300 and 500 mb) in band 3 verified highest with R? 
mameene between 0.82 and 0.91. “The previeustyenotred) frend 
in R? continues at the 500-mb level of independent data with 
R? meLeasiang poleward from 0.7/6 am Bamds) toe0. 91) dn band 
3. At 300 mb, this band distribution of Re (3) is not so 
apparent since values close to 0./2 are associated with 


bands 1 and 2 with a mean value of Ee 280082 (existing 


band 3. 


B. COMPARISONS OF FNWC INDEPENDENT VERIFICATIONS 
Comparisons of the FNWC independent tests (Table 12, 
right side) bring out primarily the same results which have 

been developed from dependent base samples. For example, 


~w 


“the FNWC overall Ro (5) ~ O.78 tay (be compared with 
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Re (3) ues” tnepeneral. the independent nq (3, indep) is 
Beeose to the values found for R= (3,dep). umeeee ple awiscn 
earlier results, band 2 has a relative minimum value of 

Ry (indep) of magnitude 0.6760. The corresponding results 
for Ro (5, indep) relative to Re (3, indep) show that the for- 
mer averages, band-for-band, between 6 and 8&4 hizher than 
the latter. Comparison of FNWC and NESC independent tests, 
both at 500 mb, shows higher mean specification of the 
latter by close to 5% in favor of NESC over FNWC independent 
Meritication in bands 2 and 3. The NESC verification sta- 
Dbality seems to be dominant in bands 2 and 3 with values of 
R? ~ 0.85. Here again, the geographic trend previously 
noted is present. 


As has been described before, Ro (5, indep) also has a 


band-2 minimum but at a level of 10Z4 higher than Ro (3,indep). 


36 





relel ye 
retlal’é, JE 


qH7C 
ole G4a bl 


aH 
reli olial 


tela ite 
qHcL 


ols OU/ 
ted ae 


YH? 
YHCT 


GWIL 
90%d 


Sou ec 
Cs Foc 


See) 
Ars y) a4 


ies Gye 
WvpesS © % 


9S eS 1 
6cc2T 2 


EGIL IL (6 
S669 ¢ 


© 160i 
COle c 


woudd 


*@Ls 


T906°0 
Sos 10 


C8E8° 0 
€678°0 


So 6.50 
T6LL°0 


aTdues Juspusdopurt 


98S56°0 
0006°0 


S773 0 
0 Gee. 


WG tO 
eS 0 


ua 
C 


B75) 9 
CUS 


OGa97 2 
EID Gee 


Orso s 
SS6E°S 


OOo etc. 
Gol Tes 


007C v7 
eG Bes 


SOS a7 = 
Sieies b- 


Oy Foe V30S 37> 
Foo 9 Chee VeH 
iG. Cece G— 
ere S Geol? 66> 
LEO 9 [7 et. Si 
Gee@ar S seas Sh 
‘Aad “uals NVANW 
CHAUAS HO 


La 
IC 


OS 
ce 


aS 
872 


02 pettTdde zoqzemtTyso Aep-aAnog (q) 


6€ 
Gis 


Et 
103 


Ge 
6 


HZis 
ATANVS 


eTdwes yuapusdopuyt of pattdde zojewtTyss aseq Aep-osetIyy (#) 


02 (OT) 


NOZ-OS 


NGG-SE 


NO7-O¢ 


NOC-OS 


NSS-SEt 


“qu OOS 328 saTdwes ysoq-juapuodoaput OSaAN 


“bap 8 OC uOpma POT Pdde euq BO SA8TNsSera LeOTast3zae32s °*8 AIAVL 


37 





38 


Os aya S704, <c BGbZ 0 t€S90°% LoLe Gs- Ly NOL-0S 
All T6SS°T Ghee 0 CES C7 BOE Ll 7S — 1G 
YH CLUES < OOS ete enue. 7 0002 0S-> OS NiGeece 
YHCT WOE 3G COa 7 0 LOO8°E Vos is— Cr 
SHC GLOe ec Bo od Geog 7 SUGGS NE vies LS NO-0Z 
EC SOMO Ge GES E70 LOSU 7 £796 29- SZ 

a—Tdues Juspusdepuy 03 pof{tTdde azoqyeuyyss aseq Aep-anog (q) 
H7.¢ ters et €L4S98°0 Gel o Osc SS 6£ NOL-0S 
piel Al Goo £868°0 Soo? 7 SCre 9S- Ge 
ely (6 CO64 °C c939 0 eS Sas Mi 9879 67- Ee Nie e=ce 
mH? OTC maa eoee © Geo 7 SOS (ne 
YH Ghou c OOnLO heeGg + 666€°74- ee qe AOe 
mic Seog. Wp US 2 ga €i7¢c 4 foe. Ly = SG 
AWTS Youur 74 tiad? dls NVAW as aqnvd 
°90uUd “dis GaAuasdo ATAIWNVS 


aTdues Juapuodapuyt of paftTdde z0j3eMTISs9N 9seq Aep-voiyL (#) 


“qu OO€ 3e saTdwes Asojy-Qusapusodspuyf, OSAN 
O62 (O01) “ba Jo voTaeoT [dae yg Jo SATNbSzr [eOFIsSTIFAS ‘*6 ATAVI 


_ 








39 


uH?7Z SASL LE 9878 °0 7T600°L GMIEG Yee i NOl-0¢ 
ane E B878°E ECL 0 866T°L TOI)  7os 9¢ 
UHV YEG Sas 96 5:2 0 GLCO OL O59 18 C= OS NigiesicG 
Oaladl LUO a ee i0) Soo 9 Cio On, oC = So 
UH?Z GLesec cLO8° 0 Coon 9 9260-70 - ES NOv-0Z 
pe gira de Lae 4c 9908 0 Cory Vil St eG 

atdues Juepusdaput o pottdde xzoyjyeutyss Aep-Anog (q) 
YHV7Z GO G4e yp ss AG) 6t27 2 CVG Ver 6€ Nos-o¢ 
ait LLG 7 OLS9°0 Gt2s 9 £89T°SE- Ge 
YH7C 6886°7C Z£T9Z°0 eee 9 WOS6 90 - Le NOC-cCE 
salve Tl SOCCER STE 10 Geo7 9 CUGL oC- Le 
YH? 2 SUS C Gao 0 ees S TES. ST ct NOv-0Z 
Ce 97 © f€66L°0 Cyac  S S96 St> AC 
AW IL Yours zu edo dis NVAW As aqnvd 
°90Ud "dis GaAYaASadO TTdWvVS 


aTdwues Juapusadsepuyz oq pottTdde azoqewy Aso Aep-svs9IyJT (e#) 


“qu OOS 38 saotdues Jsoaq-Qusapusdspuyt OMNA 
037 (OT) ‘ba fo voy eozTdde ayQ jo sqTnsot TeRoTIsTIwIS “OT ATIVI 


- = 





YHVC 


Vee aC 


S86 8) (0 


Coleg 7 


OTS8°SS- 


ayy 


40 


4HZT Tene 689° 0 9779 Y CHEZ *SS— 92 ae 
UHV 89T8°¢ LS 19, Cos0 Ss Gn OO Sia OS NiGe See 
REEL Yee eo) 7 76 fle 10) UTES Fe 6663 ° 1TS- cS 
aH7C SUSE) | S790 Geo 5G 6665 E7- LS NO¥-OZ 
UHCT L80S°¢ cSt 0 GSP8l° 7 GCo05 C7> 8¢ 

atTdues Juapuedsapuyt o3 pattdde zoRZeWMTISS 9seq Aep-ANoOg (q) 
UHV? ke Oa GZE8 0 LOC Es y 68¢L°SS~ 6€ NOL-0S 
ari 0 0678°T 79S98°0 SOGLa 7 CuSO L£S- Ce 
YH7Z SC?78 ¢ 6609°0O Las7 7 Sige 0S Le NCC-CE 
el (aa | TSSt st CULO 0 L60%7°S col?) oS We 
os Byrd CGe 7. © £66 70 LGo? 7 Cee C7- Gt NO ¥-0Z 
eUewa ie Sa. 02728 °0 Sort 7 Veg ty Ge 
AWI L Loud a 74 “AGG “aLS NVAW HZ AGS aqNnvd 
°90uUd “ais CHAXASaO ATAWVS 


atdues Juepusdapuyt o3 poyttTdde zojAyewy3ss sseq Aep-VaIY] (#) 


"qu goog 3e sotdues Ysaq-qyuapusdapuyt OMNG 
02 (OT) °bga Jo uo IeoTTdde ayqQ Jo sATnsor TeOFISTIeAIS “TT ATAVIL 





667 10 


LOGE 8 


8689°0 


70S2°0 


"BAY 
pueg 


soeL 0 


Peace 0 


G2LL°O 


SOY 
pueg 


ESE rk 10, 


See 6 918 


Sco) 


6£78°0 


€ pueg 


BEE 10 


G9 29 10 


¢T69°0 


g09550 


¢ pueg 


qu o0€ OMNA 


9€82L°0 


18h 19 


O29 c10 


¢ pueg 


76S2°0 


86SL°0 


L652 0 


¢ PpueRrg 


qu 00S OMNA 


[* (7a) QOUPTAIPA pouTeETdxd TeuoT~T_IeAJ V9Ae sSoNnTra) 
369] wofFssoeisei JQuapusdapuyt OMNA pue OSAN Fo ATeMUNS 


T9SL°0 


oo 02..0 


yp) (e 


S772 10 


T pueg 


7962°0 


7861L°0 


C7610 


T pueg 


C964) 


So cic 0 


WG Cia 


7 S22 10 


“SAV 
pueg 


C6E8°0 


ois 2 0 


6935 10 


° SAY 
pueg 


O79 G ae 


99s 0 


Sy 920 


7898°0 


g¢ pueg 


oO Ue 


vee (Ltr 218, 


Cocos 0 


SO£6°0 


€ pueg 


O10 


hoe 20 


8STL°O 


FON e «10 


79TLO 


¢ pueg 


USAN 


77Co 0 


Vove 0 


GcSO8 0 


¢ pueg 


JSAN 


84¢2°O OO€+00S 
"SAY 
Oe. 0 “SAY 
potiod 
8169°0 Aep 
-inog 
GLGL0 Aep 
-d9l1UL 

fe pueg 
6T8L°0 “BAY 
poyiog 
Z68L°0 Aep 
-inog 
941L°0 Aep 
-o90IuUl 

T pueg 


*S3[nNse2 
a ale 


41 





VI. COMPARISONS OF INDEPENDENT TO DEPENDENT RESULTS 


A summary of independent to dependent test results is 
fmven in Table 13. Three groups of tests are included. 
give first stratification (part a, Table 13) is averaged 
over all the bands and the data base modes. The outstanding 
aspect of part (a), column 2, is that it summarizes the 
comparative explained variances of Ty, and Tye The high 
explained variances for the dependent case at 500 mb may be 
noted, as well as the relatively modest decreases in R 
(between 6 and 8%) upon testing with the independent pooled 
samples. 

Pant (tb) of Table 13 averages over pressure levels and 
illustrates in summary form the distinction between the NESC 
and FNWC regression samples. These samples show the geo- 
Prapiic patterns of R? noted earlier: e.g. R (dep) has a 
tendency to increase poleward whereas Re (dep) has a band 2 
minimum. The independent test column of part (b) summarizes 
the pattern of T,.-verification by bands (but independent of 
ikewe ls)’. Upon verification, the mean-regression R? shrink- 


age of T_ averages between 4 and 6% when all bands are con- 


F 

= a ale oa celtale 
Part (c) of Table 13 Summarizes @oonewotetne Key adeas 

put forth in earlier sections.) forsexamplewwene, vartapue 


T), 05 dep) is characterized by a high value of RV (5). The 


flatter statistic decreases slightly (by about 52%) when 
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applied to independent-test data. The same remarks may be 
applied to TE ©) with slightly smaller Re (5) ~ 0.874 noted 
in the dependent case coupled with a somewhat larger shrank 
age (about 9%) after test upon independent data. The con- 
clusions regarding the comparative instability of the depen- 
dent test for TY, 03) become clear from Table 13, part (c), 
where it is seen that the values of R? increase approxi- 
mately 34 from dependent to independent tests. 

The results for TG dep) and T.(3 indep) are inconclu- 
Sive in the form presented in Table 13(c). The whole set 
Sa resm@lts Of part (¢c) indicates that little a2fsany addi— 


tional specification accuracy is obtained from the four-day 


data base in lieu of the three-day data base. 
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TABLE 13. Summary of independent to dependent test results. 


test being 
compared 


dependent 


ee 
value of R 


value of R 


(a) > 0.8957 0.8392 
T, (5) OT O7o7 0.7798 
TT, (3) 0.7298 OP Ta 35 
T, (3) ere oy 0.7201 
(b) Th band 1 0.7795 02/7546 
Tp. band 1 0.8248 Or 7 oo 
Ty band 2 0.7966 0.7701 
Th Dand) 2 0.7532 O77 7 
Th band 3 0.8620 0.8640 
Tp band 3 OS 22 O27 761 
(c) ie), 3 day 0.8989 0.8369 
T.(5), 4 day 0.8928 0.8369 
Hh ,@3); 3 day 0.7465 O) gael 
TH (3), 4 day 0.7129 0.7226 
T,(5), 3 day 0.8760 fen 5 
TE); 4 day 0.8714 Oe 7 or 
T.(3), 3 day 0.7247 0.7504 
Teo); 4 day 0.7266 0.6898 
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Vii. sCONCLUSTCNS 


ime results of —Sections) LV). eandmV lea tOrdmcom: i rman 
merdence. For each case, high specification (for Ty occurs 
at the 500-mb level. In spite of considerable shrinkage, 
there still remains a significant specification value of R 
at 300 mb. This was true also of the FNWC specifications. 

ihe reduction in R? associated with 300 mb Ty and Ty, 
specifications was maximal in band 2 where more inter-di- 
urnal variability in the temperature and associated tropo- 
pause conditions occur. The NESC specifications obtained by 
Smith et al (1970) by the 16-predictor equation also fell 
off to a relative eee at the tropopause level in their 
SIRS-A test. 

Not only is part of the observed lack of specification 
associated with the sampling of temperature extremes at or 
near the 300-mb level, but also the inter-diurnal varia- 
belity in cloud cover would affect ine mMoise contained ian 
the radiances due to the cloud-correction technique, which 
is not error-free. The errors due to cloud-correction must 
be maximal in channels 2 and 3 for which maximum dt/dln p 
weights occur in the upper troposphere. Hence, the mid- 
latitude band (band 2) which had most inter=diurnal as 
well as geographical sampling variability had the poorest 


300-mb specification. On the other hand, the 500-mb speci- 


fication and verification was adequate throughout. The 
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fact that the best 300-mb specification occurred in band 3 
(50-70N) where cloud eee are climatologically lower than 
mmethe other two bands tends to confirm this conclusion. 
Similar tests at additional levels would have provided a 
better perspective in assessing the use of the SIRS cor- 
wected radiance data in constructing vertical proriles.) «As 
meGestit. a More balanced view Jot the meutrcal resolution 
afforded by the SIRS regression techniques would have 
evolved. 

Several factors could work to improve the operational 
routine tested here: (1) to make use of a multi-level 
cloud-correction model for correction of SIRS-radiances 
(similar to that described by Smith et al 1970); (2) to 
employ a l16-predictor equation involving squares of radi- 
ances as well as linear terms; (3) to use a larger sample 
base confined to a smaller band zone, for example, 1S 
latitude overlapping bands rather than the 20° bands used 
here. 

lt is felt that with these MedifPreattons ea oh s 
operational unit attached to FNWC could profitably use the 
SIRS data in improving their update analysis over sparse 
data regions, and in determining the onset of significant 


circulation changes. 
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